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1.  Introduction 


Moving  target  indicator  (MTI)  radars  have  been  used  to  detect  and  track  targets  since  World  War 
II.  When  searching  for  targets,  a  narrow  antenna  beam  is  typically  scanned  through  a  large 
region  of  interest.  For  most  operational  radar  systems,  the  beam  is  scanned  by  rotating  an 
antenna  on  a  gimbal  (1).  For  more  recently  developed  systems,  the  beam  is  scanned 
electronically  with  a  phased  array  antenna  using  switches,  ferroelectrics,  or  other  active  devices. 
Phased  array  antennas  can  be  expensive  and  consume  large  amounts  of  power.  Another  less 
expensive  approach  for  scanning  a  radar  beam  is  to  use  a  frequency  scanning  antenna.  Systems 
using  frequency  scanning  antennas  trade  good  angular  resolution  for  poor  range  resolution  (2). 
The  goal  of  this  design  is  to  have  a  low-cost,  low-power  system  with  both  good  range  and 
angular  resolution. 

Numerous  techniques  have  been  developed  to  reduce  the  cost,  size,  and  power  requirements  of 
electronically  scanning  antennas  (ESA).  One  technique  to  reduce  the  cost  is  to  aperiodically 
place  a  smaller  number  of  antenna  elements  across  an  aperture  (5).  This  technique  has  been 
shown  to  be  feasible  based  on  performance  metrics  such  as  peak  sidelobe  level  (4).  I  am 
proposing  a  similar  concept  but  different  implementation  for  reducing  the  cost  of  the  system. 

My  proposal  is  a  radar  system  architecture  based  on  a  frequency  scanning  antenna  with  delay 
lines  whose  lengths  are  deterministic,  but  randomly  selected,  in  combination  with  digital  signal 
processing  algorithms.  This  architecture  uses  temporal  diversity  rather  than  spatial  diversity  to 
reduce  the  cost  and  improve  the  performance  on  the  system. 

Conceptually,  orthogonal  beam  patterns  at  the  receive  antenna  will  be  generated  at  each 
frequency  with  varying  levels  of  energy  spread  spatially  across  the  single  element  antenna 
pattern.  Then,  signal  processing  algorithms  based  on  techniques  such  as  least  squares  and 
maximum  likelihood  will  be  used  to  detect  and  track  the  targets  of  interest.  From  information 
theory,  the  maximum  information  content  of  a  signal  is  obtained  when  the  data  have  a  uniform 
distribution  and  minimum  information  content  is  obtained  when  the  data  have  a  delta-like 
distribution  (5).  Traditionally,  antenna  patterns  are  delta-like  functions,  which  are  scanned 
across  a  field  of  view.  More  information  about  the  targets  should  be  obtained  using  a  series  of 
patterns  that  have  distributions  more  closely  resembling  a  uniform  distribution. 

This  report  examines  the  feasibility  of  an  architecture  based  upon  a  one-dimensional  (1-D)  linear 
array  feed  with  random  length  lines.  The  major  advantages  of  the  architecture  are  a  simple  radar 
RF  design,  an  antenna  with  no  active  elements,  and  a  data  acquisition  system  comprised  of  a 
single  slow-speed  analog-to-digital  converter  (ADC).  The  major  disadvantages  of  the 
architecture  are  the  signal  processing  is  much  more  complicated,  the  effective  signal  processing 
gain  associated  with  the  antenna  may  be  reduced  relative  to  a  similarly  sized  antenna,  and  there 
are  potential  ambiguities  associated  with  target  localization. 
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2.  Radar  System 


A  radar  architecture  with  a  single  input  and  single  output  was  developed  with  the  goal  of  having 
near  orthogonal  beams  as  a  function  of  frequency  in  the  receive  antenna.  The  concept  supports 
arbitrary  waveforms  and  arbitrary  antenna  placement;  however,  to  simplify  the  signal  processing, 
a  radar  system  based  upon  a  continuous  wave  (CW)  stepped-frequency  waveform  was  analyzed. 
Figure  1  shows  a  block  diagram  of  the  radar  system.  The  radar  signal  is  transmitted  through  a 
single  antenna  and  received  with  an  array  of  antennas  that  are  combined  with  feed  lines  of 
random  length.  The  received  signal  is  amplified,  down  converted,  and  then  smoothed  with  a 
band-pass  filter  (BPF).  Only  signals  from  targets  moving  at  a  sufficient  velocity  will  pass  to  the 
ADC.  The  single  element  antenna  patterns  for  both  transmit  and  receive  antenna  elements  are 
identical. 


Figure  1.  Block  diagram  of  radar  system  with  feed  lines  for  the  receive 
antenna  whose  length  are  selected  randomly. 

Conceptually,  the  proposed  radar  system  contains  enough  information  to  track  moving  targets  in 
stationary  clutter.  The  velocity  of  the  targets  can  be  determined  by  performing  a  discrete  Fourier 
transform  (DFT)  on  the  time  samples  collected  at  each  frequency.  Conceptually,  the  range 
information  on  moving  targets  can  be  obtained  by  performing  a  DFT  across  the  Doppler- 
processed  results  after  motion  compensation  and  windowing.  Conceptually,  the  angle  of  the 
targets  can  be  determined  using  techniques  such  as  matched  filter  processing  or  least  squares 
estimation.  To  estimate  target  parameters,  stochastic  optimization  techniques  need  to  be 
developed.  This  report  studies  the  feasibility  of  estimating  the  target  parameters  by  analyzing  a 
cost  function  based  upon  the  Euclidian  norm. 
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3.  1-D  Antenna  Patterns 


Analysis  of  the  1-D  antenna  patterns  provides  insight  into  the  more  complieated  system 
identilieation  problem.  1-D  antenna  patterns  were  simulated  and  analyzed  for  a  radar  system 
similar  to  the  bloek  diagram  shown  in  figure  1 .  For  stationary  targets,  the  BPF  shown  in  figure 
1  was  removed.  The  geometry  associated  with  a  1-D  array  with  N  antenna  elements  for  the  pth 
target  is  shown  in  figure  2. 


Figure  2.  Two-dimensional  (2-D)  geometry  for  a  1-D  linear  array  for  a  stationary  target  in 
the  far-field. 


For  a  step-frequency  waveform,  the  returned  normalized  electric  field  from  P  stationary  point 
scatter  targets  in  the  x-y  plane  is  given  by 


E{m)  =  YUraA^ 


In 


fo  +  mAf 


RTp+rp(n) 


(1) 


where  E  (m)  is  the  normalized  electric  field  for  the  mth  frequency  (0  to  M-\),p  is  the  target 
number  (1  to  P),  n  is  the  antenna  element  number  (0  to  N-V),  Ap  is  the  amplitude  associated  with 
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the  pth  target, /o  is  the  earrier  frequeney,  and  A/  is  the  frequeney  step  size,  c  is  the  speed  of  light. 
Rtp  is  the  range  from  the  transmit  antenna  to  the  pth  target  and 

r  p{n)  =  Rfip  +  nLco'&O  p  + d„  + K  (2) 

is  the  range  from  the  pth  target  to  the  ADC  for  the  nth  antenna  element.  Rop  is  the  range  from  the 
left  most  reeeive  antenna  element  in  the  array  to  the  pth  target,  L  is  the  separation  between 
antenna  elements  in  a  linear  array,  and  0p  is  the  angle  to  the  pth  target  from  the  reeeive  antenna. 
Finally,  dn  is  length  of  the  nth  delay  line  and  K  is  an  arbitrary  eonstant. 

From  measurements  of  E(m),  one  ean  estimate  the  angle  of  arrival,  0p  ^  for  eaeh  target.  A 
standard  approaeh  to  estimating  0p  is  to  use  a  bank  of  matehed  filters.  The  matehed  filter  output 
ean  be  deseribed  as  an  inner  produet  between  the  measured  data  and  the  seleeted  matehed  filter 
as  shown  below 

A[o^  =  {E{m),W0j,m)^  (3) 

9 

where  Wspn)  is  the  normalized  eleetrie  field  for  a  single  point  target  at  an  angle  of  0q,  range  of 

zero  and  ealeulated  using  equation  1 .  It  is  desired  to  have  the  maximum  response  in  equation  3 
when  the  data  and  the  matehed  filter  are  at  the  same  angles  and  zero  when  they  are  at 
signifieantly  different  angles.  The  response  of  the  mateh  filter  ealeulated  for  the  pth  target 
before  the  summation  over  frequeney  is  given  by 


M-\  ❖ 

2:  E{m)w  0  (m)  = 

m=0 


where  C  is  a  eonstant  and  *  is  the  eomplex  eonjugate.  After  some  algebraie  manipulation  of 
equation  4,  we  obtain  the  following; 
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M-\  *  .^1  /O  +  wtA/'Y  . 

X  E{m)w  elm)  =CYZ-le^^'^  (^^p+^o;,) 

m=0  ^  V  y 

L(cos6’p-cos6’Jn)  + 

(L[co^ep-co^e,)n  +  {dn-dn^i)-Lco^ep)  + 
(l  (cos  6'^  -  cos  6',)  n  +  (flf „  -  J„+2)  -  2A  cos  6'^)  + 

dN-\)-iN-\)L  cos  Op) 


Ideally,  the  terms  on  the  left-hand  side  of  equation  5  should  add  to  zero  or  near  zero  for  Op  Oq. 
This  would  result  in  a  beam  pattern  with  no  or  small  side  lobes.  Now,  the  problem  is  to  seleet 
optimal  lengths  for  the  delay  lines  to  minimize  the  sidelobes  ealeulated  from  equation  5. 
Obtaining  this  goal  is  a  eomplex  nonlinear  problem.  One  approaeh  is  to  use  delay  lines  with  a 
random  statistieal  distribution,  and  then  find  solutions  that  minimize  the  expeeted  value  of  a  eost 
funetion.  Using  this  approaeh,  it  is  observed  that  the  first  term  in  equation  5  is  independent  of 
the  delay  line  length,  so  its  expeeted  value  eannot  be  set  to  zero  by  randomly  varying  the  lengths 
of  the  delay  line.  An  unbiased  solution  ean  be  obtained  by  seleeting  the  probability  distributions 
of  the  delay  lines  so  that  the  bias  from  the  remaining  terms  in  equation  5  eaneels  the  first  term. 
However,  the  realization  of  this  solution  is  not  obvious. 

Another  solution  is  to  ignore  the  bias  ereated  by  the  first  term  and  then  statistieally  minimize  the 
eontributions  from  the  remaining  terms.  A  simpler  eriterion  to  optimize  the  length  of  the  delay 
lines  is  to  minimize  the  expeeted  value  of  the  power  in  the  side  lobes  from  the  matehed  filter 
proeessing.  The  expeeted  value  of  the  eontributions  from  all  the  terms,  exeept  the  first,  will  be 
zero  if  the  probability  distributions  of  the  phase  differenees  eaused  by  the  random  delays  are 
uniformly  distributed  from  0  to  2%  for  all  eombinations.  However,  this  distribution  is  not 
realizable  in  hardware.  If  the  delay  lines  are  very  long,  the  distributions  will  be  approximately 
uniformly  distributed  due  to  the  modulation  of  the  phase  by  2%.  This  oeeurs  when  the  length  of 
the  delay  lines  satisfies 

(rf,-rf„)A//c»l  (6) 

where  n=0,. .  .N-1  and  m=0,. .  .N-1  for  all  m  and  n  exeept  when  m=n..  This  is  not  a  requirement, 
but  should  provide  a  best-ease  seenario  for  an  approaeh  based  on  statistieal  optimization.  This 
setup  results  in  sidelobe  levels  (SL)  proportional  to 
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or  _  iljl 

2^m=0  2^n=0  ^ 


fO  +  niAf 


V 


J 


i^L  {cos  0p-cos0^)n^ 


SL  =  Tt^ 


1-J^^ 


1_J2^ 


^  fO  +  mAf  ^ 

^  ^  jL{cos0p-cos0gjN 

^fO  +  mAf  ^ 

^  c  J^(cos6»^-cos^J 


(V) 


SL  =  czt:: 


m=0 


sin(;z-  /Q  +  V  (cos  6>p  -  cos  0^  N) 


sin(;r 


fO  +  mAf 


L{cos0p-cos0^)) 


and  a  signal  strength  proportional  to  NM.  This  result  is  eomplieated  and  is  diffieult  to  visualize. 
Thus,  I  ran  simulations  to  gain  more  insight  into  the  problem. 


I  used  Monte  Carlo  simulations  to  investigate  the  impaet  of  the  line  lengths  on  the  1-D  antenna 
pattern  using  two  approaehes.  The  first  approaeh  is  based  on  mateh  filter  proeessing  as 
previously  deseribed.  The  seeond  approaeh  for  estimating  the  angle  of  a  target  is  based  on 
minimizing  the  weighted  squared  error  between  a  desired  response  and  a  ealeulated  response. 
The  weighted  squared  error  in  integral  form  is  given  by 


where 


tt/I 


£=  1 


-7r/2 


/^^(0)-s^v{0)  T{0)d0 


{0)=  Z  e{-jk{nL cos 0+dn)), 

n=0 


2^(/o+fflA/) 

c 


(8) 

(9) 

(10) 


and  Bd  (0)  is  the  desired  antenna  pattern,  0  is  the  azimuth  angle,  5  is  a  veetor  of  beam  forming 
eoeffieients,  r(6)  is  the  single  element  pattern  for  a  radiator,  and  v(^is  a  M  x  1  array  manifold 
veetor  that  is  a  funetion  of  frequeney  step  number  rather  than  antenna  element  number  (d).  The 
error  is  minimized  by  breaking  the  integral  in  equation  8  into  a  summation  over  /  intervals,  and 
then  approximating  eaeh  interval  as  a  single  digital  sample.  Now,  the  solution  to  minimizing  the 
weighted  squared  error  is  given  by 

^  =  Wv  (11) 

9 

where 

H  =  [v{0,)  v{0,)  ...  v(^0]’  (12) 


and  IF  is  a  diagonal  matrix  with  elements  given  by  r(6’/)  (7). 
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4.  Simulation 


I  performed  simulations  to  investigate  the  1-D  antenna  patterns  of  stationary  targets  and  the  cost 
function  of  moving  targets.  I  selected  the  parameters  for  the  simulation  based  on  the  competing 
goals  of  having  a  low-cost  system  and  near  orthogonal  antenna  patterns.  Simulations  were 
performed  with  a  carrier  frequency  of  10  GHz,  a  bandwidth  of  1  GHz,  64  frequency  steps,  64 
antenna  elements  with  A,/2  antenna  element  spacing  and  Gaussian  beam  patterns  centered  at 
broadside  with  one-way  3-dB  beamwidths  of  120°,  and  delay  lines  lengths  that  were  randomly 
selected  with  a  uniform  distribution  between  0  and  20  X.  Figure  3  shows  the  antenna  patterns  for 
the  first  frequency  at  10  GHz,  the  second  frequency  which  is  incremented  by  15.6  MHz,  and  the 
fifth  frequency,  which  is  incremented  by  4  x  15.6  MHz. 


Figure  3.  The  1-D  receive  antenna  patterns  for  the  first,  second,  and  fifth  frequency  step. 

The  orthogonality  of  the  antenna  patterns  can  be  analyzed  by  performing  a  single  value 
decomposition  (SVD)  on  the  observation  matrix  calculated  in  equation  12.  Figure  4  show  the 
eigenvalues  associated  with  the  simulation  parameters.  The  results  indicate  that  eigenvalues  are 
near  zero  around  the  30*  eigenvalue.  Clearly,  the  antenna  patterns  are  not  orthogonal.  The  rank 
of  the  matrix  can  be  improved  by  using  more  elements,  longer  delay  lines,  and  more  bandwidth, 
but  it  is  never  full  ranked  for  a  realizable  system. 
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Figure  4.  Eigenvalues  of  the  observation  matrix  ereated  from  the  different  antenna  patterns. 

One-dimensional  antenna  patterns  were  simulated  using  algorithms  based  on  matched  filter 
processing  and  weighted  least  squares  (WSLQ).  The  matrices  required  for  the  WLSQ  fit  were 
generated  by  sampling  the  array  manifold  vector  in  angle  every  7i/512  radians  and  the  matrix 
inverse  was  computed  in  MATLAB  using  the  pseudo  inverse  function.  The  results  are  shown  in 
figure  5.  The  red  curve  shows  the  results  for  a  WSLQ  fit  to  the  desired  antenna  pattern  shown  by 
the  blue  curve,  and  the  black  curve  is  generated  using  conventional  or  a  matched  filter  approach. 
As  expected,  the  sidelobes  are  high. 


Figure  5.  Simulated  1-D  antenna  pattern  using  matched  fdter  and  WLSQ  approaches. 
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The  sidelobes  levels  in  the  antenna  patterns  can  be  reduced  using  adaptive  digital  signal 
processing  techniques  that  do  not  distort  the  desired  signal  response  while  minimizing  the  signals 
from  another  direction.  Mathematically,  these  conditions  are  described  using 


w^V  (6'o)  -  0 


(13) 


w"v[e,) 

~V{04 


(14) 


where  6a  is  the  desired  angle,  Oq  is  the  desired  null  angle,  and  |  |  is  the  Euclidian  norm.  The 

corresponding  filter  coefficients  can  be  calculated  by  projecting  the  unconstrained  solution  or  the 
matched  filter  solution  on  to  the  subspace  constrained  by  equation  13  [6  -  van  trees].  The 
corresponding  projection  matrix  is  given  by 

=  V (^o)[F(^o)^F(^o)]  V{6oT  (15) 


and  the  resulting  coefficients  are  given  by 

wf^V{0j)(l„-Pgy  (16) 

These  equations  will  reduce  the  sidelobe  level  at  a  desired  angle.  To  reduce  all  the  significant 
sidelobes,  the  null  angle  can  be  scanned  while  keeping  the  desired  angle  fixed,  and  then  finding 
the  minimum  response  as  shown  below: 

Y  {0)=mm\w^jV{e],  y  =  !.../.  (17) 

where  Wj  are  the  coefficients  calculated  for  the  jth  null  angle  and  v(^  is  a  M  x  1  array  manifold 
vector  for  a  given  angle  that  is  a  function  of  frequency  step  number  (5).  This  procedure  requires 
that  only  one  signal  be  present.  For  applications  where  the  targets  are  sparse  and  the  signals  can 
be  separated  in  Doppler,  this  may  be  a  reasonable  assumption.  If  there  is  more  than  one  target, 
then  the  additional  targets  need  to  be  located,  and  then  nulled,  before  the  new  target  is  detected. 
This  process  requires  an  additional  degree  of  freedom  for  each  additional  target  and  slightly 
degrades  the  performance  of  the  algorithm.  The  results  from  figure  5  were  reprocessed  using 
equation  17. 

One-dimensional  antenna  patterns  were  simulated  for  signal-to-noise  ratios  (SNRs)  of  10,  20, 
and  25  dB,  as  shown  in  figures  6  and  7.  The  SNRs  were  calculated  for  a  single  antenna  element. 
The  array  of  null  angles  went  from  -nil  to  nil  in  7t:/256  increments,  with  the  region  at  the  desired 
angle  ±  7r/128  excluded.  The  antenna  patterns  were  calculated  at  angles  going  from  -nil  to  tiH 
in  7t:/256  increments,  but  were  offset  by  7t:/512. 


9 


Figure  6.  The  1-D  antenna  patterns  simulated  for  SNRs  of  10,  20,  and  25  dB. 


Figure  7.  Zoomed  in  1-D  antenna  patterns  simulated  for  SNRs  of  10,  20,  and  25  dB. 


The  results  indieate  that  the  maximum  peak  sidelobe  level  is  ~15  dB  below  the  signal  of  interest 
for  an  SNR  of  10  dB.  The  performance  improves  with  an  increasing  SNR.  The  Rayleigh 
criterion  for  angular  resolution  for  a  linear  array  at  broadside  is  given  by 


A(9 


f  \ 

c 


[L{f,+BI2)j 


(18). 


For  the  parameters  selected,  this  corresponds  to  an  angular  resolution  of  ~0.03  radians  or  1.7°. 
This  angular  resolution  was  almost  achieved  in  the  plots  shown  in  figure  7. 

The  simulation  results  indicate  that  reasonable  1-D  patterns  can  be  generated  for  stationary 
targets.  The  sidelobe  levels  in  the  1-D  patterns  improve  as  the  number  of  elements  and  the 
fractional  bandwidth  is  increased.  The  performance  also  improved  as  the  number  of  frequencies 
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is  increased,  but  levels  off  as  it  approaehes  the  number  of  antenna  elements.  For  the  parameters 
seleeted,  reasonable  sidelobe  levels  required  a  relatively  high  SNR. 


5.  Cost  Function 


A  radar  system  based  upon  this  arehiteeture  needs  to  deteet,  and  then  estimate  the  loeation  and 
veloeity  of  targets  in  its  field  of  view.  The  relationship  between  the  parameters  required  for 
traeking  and  the  target  returns  is  eomplex  and  eannot  be  estimated  using  eonventional  radar 
signal  proeessing  algorithm.  The  estimation  problem  ean  be  formulated  as  a  nonlinear 
optimization  problem  of  a  eost  funetion.  The  eost  funetion  ean  be  minimized  using  a  variety  of 
teehniques  depending  upon  the  applieation.  Genetie  algorithms,  simulated  annealing,  or  partiele 
swarm  algorithms  ean  be  used  to  minimize  eost  functions  with  multiple  loeal  minimums  (P). 
These  teehniques  ean  be  eombined  with  traditional  approaehes  sueh  as  gradient  deseent, 
Newton’s  method,  or  the  Nelder-Mead  method  onee  the  global  minimum  is  loeated  {10). 


One  possible  eost  funetion  to  minimize  the  Euelidian  or  Li-norm  is 


min(Ap,Xp,Vp,^p) 


D(m,k)-^(  D(m,k,Ap,Xp,Vp,^p) 


p=\ 


(19) 


where  Ap  is  the  amplitude  of  the  pth  seattering  eenter,  Xp  is  the  position  of  the  pth  scattering 
eenter,  Vp  is  the  velocity  of  the  pth  seattering  eenter,  ^p  is  the  phase  of  the  ith  seattering  eenter, 
P  is  the  total  number  of  target  scattering  centers,  D(m,k)  is  the  DFT  of  n  data  values  measured  at 
the  m'^  frequeney  step,  k  is  the  Doppler  bin  number,  and  D(m,k,Ap,Xp,Vp,^p)  is  the  predieted 
value  of  D(m,k).  This  is  a  eomplieated  optimization  problem,  but  for  some  applieations,  the 
number  of  targets  in  each  Doppler  bin  should  be  zero  or  one,  and  for  many  applieations  there  is 
only  one  or  a  few  seattering  eenters  per  target.  This  should  simplify  the  problem  eonsiderably. 

Two-dimensional  shoes  of  the  eost  funetion  desoribed  in  equation  19  were  generated  for  a  target 
with  an  initial  position  of  (5,  0)  m  and  a  veloeity  of  (15,  15)  m/s  in  the  x-y  plane,  with  no  phase 
offset,  acceleration,  or  noise  using  the  same  radar  parameters  as  previously  desoribed  for  the  1-D 
antenna  simulation.  Two-dimensional  shoes  of  the  eost  funotions  are  plotted  in  figures  8  through 
1 1  as  a  funetion  of  position  and  veloeity,  while  the  other  parameters  are  fixed  at  their  oorreot 
values.  The  plot  increments  are  at  0.8  of  the  theoretical  range  and  veloeity  resolution  derived 
from  the  Rayleigh  oriterion. 


The  theoretioal  angular  resolution  for  a  stationary  target  given  by  the  Rayleigh  oriterion  is  by 


A0 


r  c  1 


^Lcos(a)(/o+R/2)  J 


(20) 
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(21) 


where  a=atan(0/5)=l.  The  figures  are  plotted  at  increments  given  by 


Ax 


0.8 


^  c  ^ 


2B 


Ay  =  O.S{rA0) 


(22) 


Ar=o.8(^) 


(23) 


where  c  is  the  speed  of  light,  B  is  the  bandwidth  of  the  transmitted  signal,  r  is  the  range  (5  m), 
and  T  is  the  time  for  each  step-frequency  waveform. 


Figure  8.  Cost  function  for  a  single  target  located  at  (5,0)  m. 
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Figure  9.  Cost  function  zoomed  in  for  a  single  target  located  at  (5,0)  m. 


Figure  10.  Cost  function  for  a  single  target  with  a  velocity  vector  of  (15,15)  m/s. 
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Figure  11.  Cost  function  zoomed  in  for  a  single  target  with  a  velocity  vector  of  (15,15)  m/s. 

For  these  simulated  parameters,  the  target  is  clearly  localized  in  velocity  and  downrange  with 
resolution  equal  to  the  Rayleigh  criteria,  but  the  minimum  associated  with  the  cross-range 
position  of  the  target  has  a  relatively  low  separation  between  competing  local  minimums.  The 
localization  of  the  cross-range  position  will  be  more  sensitive  to  noise,  but  the  global  minimum 
is  centered  around  the  correct  value  and  is  slightly  wider  than  the  value  predicted  by  the 
Rayleigh  criteria.  This  result  is  expected,  since  the  stationary  1-D  antenna  patterns  also  had 
large  sidelobes.  The  sidelobes  of  the  cost  function  plots  can  be  improved  by  changing  the  radar 
and  waveform  parameters.  Increases  in  the  percent  bandwidth  of  the  waveform,  the  length  of  the 
time  delays,  and  the  number  of  elements  will  reduce  the  sidelobes  in  the  plots  of  the  cost 
functions.  Also,  performance  improvements  can  be  achieved  by  extended  this  concept  from  a 
single  array  with  one  ADC  to  multiple  subarrays  with  a  dedicated  ADC  for  each  subarray.  These 
results  indicate  that  it  is  possible  to  estimate  the  position  and  velocity  of  a  point  scatter  target 
from  2-D  slices  of  the  proposed  cost  function  with  resolution  close  to  the  Rayleigh  criterion. 


6.  Conclusions 


A  simple  radar  architecture  with  a  single  input  and  single  output  was  developed  that  supports 
three-dimensional  (3-D)  electronically  scanning  of  moving  targets  with  high  spatial  resolution. 
The  architecture  is  based  on  a  frequency  scanning  antenna  with  random  length  feed  lines.  In  this 
study,  I  investigated  the  feasibility  of  the  architecture  to  perform  target  localization  in  position 
and  velocity  by  analyzing  1-D  antenna  patterns  and  a  cost  function  based  upon  an  L2-norm.  A 
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single  realization  of  the  system  was  selected,  simulated,  and  analyzed.  The  results  indicate  that 
the  concept  is  feasible,  but  requires  a  high  SNR  and  complex  signal  processing. 

For  the  selected  parameters,  the  1-D  antenna  patterns  generated  using  matched  filter  and 
weighted  least  squares  processing  had  high  sidelobe  levels.  To  improve  the  results,  the  sidelobes 
were  nulled  over  all  angles  except  the  desired  angular  region.  A  final  antenna  pattern  was 
formed  by  taking  the  minimum  value  of  each  antenna  pattern  at  every  angle.  This  resulted  in 
antenna  patterns  with  a  resolution  predicted  by  classical  array  theory,  but  it  required  increased 
signal  processing  and  was  susceptible  to  interference  by  other  targets.  Two-dimensional  patterns 
can  be  generated  using  the  same  techniques. 

Target  localization  will  require  nonlinear  optimization  of  a  cost  function.  Two-dimensional 
slices  of  a  cost  function  based  upon  an  L2  norm  indicated  that  global  minimum  exists  for 
localizing  simple  targets.  However,  the  y-position  estimate  had  a  relatively  small  separation 
between  other  competing  local  minimums.  The  range  and  velocity  resolution  seen  in  the  plots  of 
the  cost  function  were  consistent  with  the  range  and  velocity  resolution  calculated  from  the 
Rayleigh  criterion. 

For  the  system  to  be  practical,  the  solution  to  the  localization  algorithm  should  be  calculated 
from  the  radar  platform.  Doppler  detections  that  exceed  a  given  threshold  could  be  wirelessly 
transmitted  to  a  central  computer  for  more  processing.  This  process  would  reduce  the  cost,  size, 
and  power  requirements  of  the  radar  front  end.  The  power,  size,  and  cost  of  the  central  computer 
system  would  be  larger  than  traditional  radar  systems,  but  it  could  serve  multiple  roles  and  be 
located  in  a  secure  location  with  adequate  resources. 
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List  of  Symbols,  Abbreviations,  and  Acronyms 


1-D 

one-dimensional 

2-D 

two-dimensional 

3-D 

three-dimensional 

ADC 

analog-to-digital  converter 

BPF 

band-pass  filter 

CW 

continuous  wave 

DFT 

discrete  Fourier  transform 

ESA 

electronically  scanning  antenna 

MTI 

moving  target  indicator 

SL 

sidelobe  level 

SNR 

signal-to-noise  ratios 

SVD 

single  value  decomposition 

WSLQ 

weighted  least  squares 
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